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Gluodynamics in external field in dual superconductor approach 
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Abstract 

We show that gluodynamics in an external Abelian electromagnetic field should pos- 
sess a deconfining phase transition at zero temperature. Our analytical estimation of the 
critical external field is based on the dual superconductor picture which is formulated in 
the Euclidean space suitable for lattice calculations. A dual superconductor model cor- 
responding to the SU{2) gluodynamics possesses confinement and deconfinement phases 
below and, respectively, above the critical field. A dual superconductor model for the 
SU{3) gauge theory predicts a rich phase structure containing confinement, asymmet- 
ric confinement and deconfinement phases. The quark bound states in these phases are 
analyzed. Inside the baryon the strings are Y-shaped as predicted by the dual supercon- 
ductor picture. This shape is geometrically asymmetric in the asymmetric confinement 
phases. The results of the paper can be used to check the dual superconductor mechanism 
in gluodynamics. 

1. At present there are two popular approaches to the problem of color confinement in 
gluodynamics. They are based on the Abelian monopole |1| and on the center vortex |0 
pictures of the gluodynamics vacuum. In this paper we discuss the Abelian monopole approach 
which suggests that confining degrees of freedom of the vacuum in an Abelian projection 
can be described as a dual superconductor. The key element of this picture is the monopole 
condensate which squeezes a chromoelectric flux to a confining string due to the Meissner effect. 
The string is an analog of the Abrikosov vortex in an ordinary superconductor while the 
Abelian monopoles are playing the role of the Cooper pairs. This picture has been confirmed 
in many numerical simulations on the lattice (for a review see, e.g., Ref. 

Here we investigate the properties of the SU{2) and SU{3) gluodynamics in the external 
electromagnetic field using the dual superconductor approach. Our study is motivated by 
the fact that the response of the vacuum of a gauge theory on external fields may provide 
interesting information about the vacuum structure. The external fields were used to study 
nonperturbative properties of QCD 0, the baryogenesis in the electroweak theory [0, features 
of various topological defects in three dimensional models 0, etc. 

A common feature of known superconductors is the Meissner effect: the superconductors 
expel relatively weak external magnetic flux from their interior. Strong enough magnetic field, 
^ext ^ j^^^ destroys the superconductivity and the superconductor goes in the normal (metal) 
state. In this paper we estimate analytically the critical electromagnetic fields which break the 
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dual superconductivity in the SU (2) and SU{'i) gluodynamics. Consequently, the confinement 
is (partially, in the case of S't/(3)) lost at these critical fields. 

We consider the external electromagnetic fields in a particular Abelian gauge, which is 
used to define the Abelian degrees of freedom. Strictly speaking, these purely Abelian fields 
are unlikely to be realized in NatureQ. Nevertheless, the response of the gluodynamics media 
on these purely Abelian fields can be used for further checks of the dual superconductivity 
hypothesis in numerical simulations of lattice gluodynamics. Since the dual superconductor 
picture describes various non-perturbative phenomena such numerical test is physically 
motivated. 

We consider the gluodynamics at zero temperature because in this case the relevant cou- 
plings of the corresponding dual superconductor models are already known. Indeed, in the 
absence of the external field the gluodynamics experiences the deconfining phase transition 
at sufficiently high temperature. On the other hand, at the critical temperature the dual 
superconductivity was demonstrated to be destroyed Thus the couplings of the dual super- 
conductor (at least, the value of the monopole condensate) must depend on the temperature. 
This dependence is not known at the time being. 

Apart from knowledge of the dual couplings, another simplification comes from the Eu- 
clidean formulation. The action of the four dimensional dual superconductor is just a trivial 
dimensional generalization of the (Helmholtz) free energy functional of the ordinary three di- 
mensional superconductor. Thus, from the point of view of static effects - such as a response 
of the superconductor on a static external field - the Euclidean dual superconductor model of 
the gluodynamics describes just an infinitely large four dimensional superconducting material. 
The electric and magnetic components of the electromagnetic field differ only by the orientation 
of the field strength tensor in the coordinate space with respect to the time axis. However, 
in the Euclidean formulation at zero temperature no distinguished time-direction exists. In 
this particular respect there is no difference between external the static electric and the static 
magnetic fields (this in no more correct in the presence of the external sources such as heavy 
quarks or monopoles). Therefore in this article we are using the terminology "electromagnetic 
(EM) field" for this particular case. 

Yet another simplification is due to the fact that the vacuum of SU (2) and SU{'i) Yang- 
Mills theories in the Abelian projection is known to be close to the border between type-I and 
type-II dual superconductors |T0[. At the borderline - called also "the Bogomol'ny limit" ||1T 



- analytical results for the string tension are available. Below we consider the dual models for 
both theories in the Bogomol'ny limit. We assume that the external field does not change the 
couplings of the dual superconductor modelQ. 

2. Let us first consider the SU{2) gauge theory in the 4D Euclidean space. The infrared 
properties of the vacuum of this model can be described by the Abelian Higgs (or, Ginzburg- 
Landau) Lagrangian: 

^gl[B, $] = ^F^l + l\D,{B) $1^ + a(|<|.|^ - 77^)' , (1) 



^Note, however, that a general non-Abehan field may have a non-zero projection on the Abelian subspace 
discussed in the paper. 

^This assumption works well in the macroscopic (Ginzburg-Landau) description of the ordinary supercon- 
ductivity. 
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where F, 



fii/ 



diyBfj, is the field strength for the dual gauge field 5^, $ is the monopole 
field with the magnetic charge gM and = + iguB^ is the covariant derivative. The 
gauge field B^ is dual to the third component of the gluon field in an Abelian gauge. The 



model possesses the dual U{1) gauge symmetry, B^ ^ B^ — (9^a, $ — >• e^'J^"- $. The form of 
the potential implies the existence of the monopole condensate, |($)| = Tj with rf > 0, and, 
consequently, non-zero masses of the dual gauge, niB = gr], and monopole, m$ = 2\/2A?7, fields. 

The Bogomol'ny limit corresponds to a region of the coupling space where the masses of the 
monopole and gauge fields are the same. In our notations (|l]) the Bogomol'ny limit is defined 
by the condition 



(2) 



The properties of the Ginzburg-Landau model ([^) are very well known. Below we briefiy 
derive the value of the critical field in the U{1) model (|l]) and later we apply this method to a 
more complicated [[/(l)]^ case corresponding to the SU{3) gluodynamics. 

Consider the four-dimensional sample of the (dual) superconductor occupying half-space, 
X2 ^ 0. Let us apply the constant external EM field F^^* = e^v^i^H'^^^ to the boundary of the 
superconductor. A weak external field penetrates inside the sample up to the distance ~ m 



B ■ 



The screening of the field (the Meissner effect) is reahzed due to the induced superconducting 
current. 



= 53m($*Z}^(fi)$) 



d^ip + gAiBf, 



(3) 



where we have set $ = l^le**^. The current is parallel to the boundary of the superconductor. 



The monopole kinetic term in Eq. (|T]) can be written as 



i$iv 



Clearly, 



a non-zero current provides an additional positively-defined term in the Lagrangian (oc |$P). 
As a result, the external field lowers the value of the monopole condensate. We are looking for 
the critical value of the EM field which destroys the monopole condensate and, consequently, 
confinement. 

We disregard the quantum fiuctuations of the fields in the model (|1]) treating this system 
classically. To derive the critical EM field it is convenient to rewrite the action of the model ([1|) 
as a two-dimensional integral. The first direction of the two-dimensional plane is obviously 
the depth of the dual superconductor, X2, while the second is given by the direction of the 
current (|^). We choose oc i and rewrite the model (Q) in the (1,2) plane. The first term 
of eq. (|I|) is = while the second term can be rewritten with the help of relations: 



E 



=1,2 



(4) 



The last term in this equation can also be represented as i^*[Di, D2]^ 
the following action in the (1, 2) plane: 



S, 



GL 



L3L4 



2 2 



V 



guH |$| . We get 
+ Sj , (5) 



2i 
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where Li is the (infinite large) length of the dual superconductor in ith direction and Sj is the 
action of the surface current: 



Sj = L3L4 / d^Xeaf3daJf3 = L1L3L4 ■ Ji{x2 = 0) . 



(6) 



Using condition 
1 



we rewrite Eq. further: 



S 



+ 



T Sfiux + Sj . 



(7) 



The sign in this equation is dictated by the "fiux" action, Sfiux = gMLzL^rf j d?xH/2, which 
must be positive. We choose H > and get two Bogomol'ny equations which minimize the 
action (|^: 

The second equation gives the value of the monopole condensate at the boundary 



mx2 = o)\' 



V 



M 



The condensate disappears at the critical value of the EM field, H"^^^ = Her = dMV^/'^- At this 
value the superconducting current at the boundary vanishes, Ji = 0, giving Sj = 0. The total 
action is given by the fiux contribution, S = Sfiux = Vol ■ gj^rj^/i which is the sum of the free 
energy of the normal state, Sn = Sgl[B = 0, $ = 0] and the free energy of the EM field, Sh- 
Sn = SH = Vol ■ glr^ys. 

In the Bogomol'ny limit the tension of the string spanned on trajectories of the fundamental 
charges can be evaluated exactly^ [|Tl|, |12 



(gM) and electric (^f) charges, gMg 
of the string tension: 



a = nrj . Using the Dirac relation between magnetic 
27r, we get the exact value of the critical EM field in terms 



gHcr/a = 1 , for SU(2) . 



(9) 



This equation can be used in numerical simulations to check independently the closeness of the 
gluodynamics vacuum to the boundary between type-I and type-II superconductors. 

The nature of the phase transition at the critical EM field can be understood as follows. 
The EM field, applied to the boundary of the media with the SU (2) gauge fields, lowers the 
value of the monopole condensate according to Eq. (^. For the external fields H'^'^^ < Her 
the monopole condensate is non-zero, |$(x2 = 0)| > 0. This implies, that the dual photon 
mass at the boundary is non-zero as well, mB{x2 = 0) = gM\^{x2 = 0)|. Due to the Meissner 
effect the external field diminishes as we go deeper into the media. At the distances of the 



order of the correlation length. 



the monopole condensate restores at bulk value, $ ~ 



Consequently, at these distances the dual photon mass is restored as well, ~ gMV and the 



external EM field diminishes exponentially, H{x2) 
is present in the bulk of the system. 



~ e-'^BXi fo^ > rpj^^^g^ no EM field 



^We stress that this result is obtained at zero temperature in the absence of the external fields. 
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Now suppose that the external field reaches its critical value, H^^^ = Her- Then both the 
monopole condensate and the dual photon mass are zero at the boundary. The last fact implies 
that the field penetrates inside the media for an infinitesimally small distance 6x2 without any 
suppression because the Meissner effect is absent. At the distance 6x2 the situation repeats 
again: the condensate and the photon mass are zero and the field penetrates deeper into the 
media for another infinitesimal step, etc. Thus, the external field "eats" the condensate step 
by step and finally the condensate disappears in the whole space. Due to this mechanism the 
EM field originated at the boundary can destroy the monopole condensate (and, consequently, 
the confinement) in the bulk. 

3. Now let us consider the Lagrangian of [?7(1)]^ Higgs model corresponding to SU{3) 
gluodynamics [ll3| : 



C = \f^^F-'^-^ + J2 [l\Di:^^^f + a(|«I>.|' - v'Y] , (10) 



1 „., ^ri,„/,\, ,2 ,9 oX^i 

i=l 



where = d^B^ — d„B^ is the field strength for the gauge fields 5°, a = 3,8, Dy 
dfi. + WM^iB^ is the covariant derivative acting on the monopole fields $i, i = 1,2,3. The e's 

are the root vectors of the group SU{3y. el = (1, 0) , ea = (-1/2, -^3/2) , eg = (-1/2, ^3/2). 
No summation over the Latin index i is implied. 

The gauge fields B^'^ are dual to the diagonal components a = 3,8 of the gluon field 
y4^. Lagrangian ([TOD respects the dual [t/(l)]^ gauge invariance: B^ ^ B'^ + 6i 
Gi + fi'Af(£^itt^ + ^fo^*)) a = 3,8, i = 1,2,3, where and are the parameters of the gauge 
transformation. The phases of the monopole fields satisfy the relation 



arg$, = 0, (11) 

' '1=1 

which plays an important role in the formation of the quark bound states within the dual 



superconductor formalism |]T^ 115 



The Bogomol'ny limit is defined by condition [^6|, |1^ 

gl/X = 16/3, (12) 
and the equations of motion in this limit are given by 

(d«±zD«)x. = 0, H^'^t^-^{\x.\'-v')=0; ^ = l,2,3, (13) 

where 

ifW = J2 ^H"", (14) 

a=3,8 

are the EM fields projected on the (3, 8)-charges of the monopole fields. One can imagine these 
three fields as the dual red, dual blue and dual green EM fields. 
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The second equation in ( ]TB| ) gives the same critical value for all components of the if*^*)- 
fields: 



i/«=if„ = 3W/4. (15) 

The critical values are equivalent due to the Weyl symmetry of the dual model |]T^, TE\ which 
states that the [[/(l)]^ Lagrangian (|iy) is invariant under the transformations of the dual gauge 
fields B^'^ corresponding to the mutual permutations of the if'-*-' fields ([l^). 



The string tension spanned between the fundamental charges (quarks) in the Bogomol'ny 
limit ([l2|) of the SU{3) gluodynamics is ||16|, [T^ a = 2'Krf. Using the Dirac quantization 



condition we get the critical field Her in units of the string tension: 

gHer/(T = 3/A, forSU(3). (16) 

When the strength of the EM component if*-*^ reaches the Her value then the expectation value 
of the corresponding component of the monopole field, $j, gets vanished. Note, however, that 
the fields if*^*^ play an auxiliary role because they are not independent according to Eq. ([I^) . 
Expressing the auxiliary fields if'-*^ in terms of the components the EM field, if'^'^, and using 
Eqs. ([T5|JT6|) we get the phase diagram depicted in Figure |I|. 



The phase diagram contains confinement (C), deconfinement {D) and the asymmetric con- 
finement phases {A). The position of the phase transition depends not only on the absolute 
value of the EM field but also on the ("color") orientation of this field in the Cartan subgroup. 
At low values of the field the model is always confining regardless of the color orientation. 
However, as the absolute value of the field is increased, the model enters - depending on color 
orientation - one of six (A12, A13, ^23, Ai, A2 or A^) asymmetric confinement phases. In the Aij 
phase the ith and jth components of the monopole field are condensed while the expectation 
value of the third component is zero. In the phase Ai the ith component is condensed while the 
others two components are not. With the further increase of the field the model either enters 
the deconfinement phase, or stays in one of the three asymmetric confinement phases, Ai, 
A2 or A3. 

4. The quark bound states in the Abelian projection approach are classified with the help 
of the states of the strings spanned between the constituent quarks. The string configuration 
in the baryon was extensively studied both in the Abelian projection approach [[1^, p!5| as well 



as in a gauge independent formalism |T^. The dual superconductor model predicts [jl5| the 
existence of the F-shaped string configuration in agreement with most of Ref. |19[ . Here we 
discuss the quark bound states in the presence of the external Abelian field. 

The confining properties of both the confinement and the deconfinement phases are stan- 
dard. In the confinement phase all three monopole fields, $j, i = 1,2,3, are condensed and 



the [t/(l)]^ model possesses three types of vortex solutions [T3, IT3, EDI which confine all quarks 



into bound states. Each vortex solution is characterized by the winding number n^, $i oc e*"*''^ 
where 0j is the azimuth angle. The winding numbers of the strings are subjected to the con- 
straint Hi + n2 + = coming from Eq. ([TT|). In the deconfinement phase all monopole fields 
are not condensed and the bound stated are not formedQ. 

^In this paper we disregard weakly bounded states which might appear due to the exchange of the pertur- 
bative gluons between quarks. We also disregard the role of the Jacobian arising in the string model. 
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I 



Figure 1: The zero temperature phase diagram of the dual Higgs model in the Bo- 

gomol'ny limit in the presence of the external electromagnetic field. The confinement and 
deconfinement phases are denoted as C and D. The asymmetric confinement phases are de- 
noted as Ai and Ajk, where subscripts indicate the components of the monopole field $j which 
are condensed in these phases. 

Let us discuss the bound states in the asymmetric confinement phases. Consider, for exam- 
ple, Ai2 phase. In this case the model possesses only two condensates, ($1,2) 7^ 0. Consequently, 
on the classical level only two types of the chromoelectric strings can be formed with ni^2 7^ 0. 
Despite the third field has zero expectation value, = 0, its phase may fluctuate and the 
(classically tensionless) 71,3 7^ string-like configurations may appear. This implies that in 
addition to the purely classical n = (1,-1,0) string configuration there exist also (1,0,-1) 
and (0, 1, —1) configurations. These additional configurations are composed from the classical 
strings with either ni 7^ or n2 7^ and a string-like (tensionless) quantum excitation with 
ns 7^ 0. The stability of the quantum components of such configurations is guaranteed by 
Eq. (O).^ 

The n = (1,-1,0) string configuration must have bigger string tension than (—1, 0, 1) and 
(0, 1, — 1) configurations. Thus in the A12 phase we expect the existence of the relatively heavy 
meson composed from red quarks and two lighter mesons made of blue and green quarksQ 
The observed color asymmetry in the meson states is caused by the breaking of the Weyl 

^Here we adopt the classification of Ref. assuming that the quarks in a meson state are connected with 
each other by a pair of the strings with winding numbers n — (1,-1,0), (—1, 0, 1) and (0, 1, —1) for RR, BB 
and GG states respectively. 
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symmetry |p!7| , ^ of the Lagrangian (|TU|) by the external field. 

The baryon state also exists in the A12 phase but it should be lighter than the baryon in the 
confinement phase. The quarks in the baryon are connected to each other by all three types 
of the string configurations. It is known [|15|] that in the absence of the external field these 
strings form a symmetric ^-shaped profile. In the A12 phase the tension of the (1, —1, 0) string 
configuration is heavier than that of (—1,0, 1) and (0, 1, —1) configurations. Therefore in this 
phase the strings in the baryon state must form an asymmetric ^-profile, see Figure |[ 
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Figure 2: The meson and baryon bound states in confinement, deconfinement and A12, Ai 
asymmetric confinement phases. The solid lines correspond to the strings existing on the 
classical level while the dashed lines represent the quantum string excitations. 

The properties of the A^ phases are different from those of the Aij phases. In particular, in 
the Ai phase the string with the winding number rii 7^ has a non-zero string tension while 
the strings with 77,2,3 7^ are classically tensionless. Thus in this phase only (1,-1,0) and 
(1,0,-1) string configurations have a non-zero string tension. This implies the existence of the 
mesons composed of red and blue quarks while the meson made of green quarks is absent. The 
existing mesons must be lighter than those in the confinement phase. All possible meson and 
baryon states in confinement, deconfinement and asymmetric confinement phases are depicted 
in Figure 

5. Summarizing, we have explored the phase structure of the SU(2) and SUi?)) gauge 
theories in the external electromagnetic fields at zero temperature. Both theories are considered 
in the dual superconductor formalism formulated in the Euclidean space (suitable for lattice 
calculations) and in the Bogomol'ny limit (as confirmed by various lattice calculations). 

We have found that the phase diagram of the dual superconductor model corresponding 
to the SU{2) gluodynamics contains confinement and deconfinement phases which are located 
below and, respectively, above the critical field. The critical electromagnetic field is analytically 
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estimated in terms of the string tension in Eq. (^). 

The phase diagram of the dual superconductor corresponding to the SU{3) gluody- 

namics contains 8 phases (confinement, deconfinement and 6 asymmetric confinement phases). 
The reach phase structure - shown in Figure |I| - comes from the dependence of the monopole 
condensate on the color orientation of the external field. This finding is supported by the obser- 
vation of Ref. that the properties of the Abelian monopoles depend on the color orientation 
of the monopole. Three of six asymmetric confinement phases {A12, A13 and ^423) contain one 
baryon and three meson states. Two of the meson states are lighter then the third one. The 
strings in the baryon form the asymmetric F-shaped profile. These phases can be still regarded 
as confinement phases since the quarks of all three colors are confined. 

The other three asymmetric confinement phases [Ai, A2 and A3) may contain only two 
light meson states while the baryon state is absent at all. The quarks carrying a particular 
(phase-dependent) color are not confined in these phases. 

Note, that our results were obtained for the Abelian external fields which are applied in a 
fixed Abelian projection^ In particular, this means that the recent results of Refs. [^| for the 
phase diagram in the external Abelian fields can not be compared with our predictions because 
these results were obtained without the gauge fixing. 

It would be interesting to check the predictions of this paper by numerical simulations 
performed in the Maximal Abelian projection of the SU{2) and SU{3) gluodynamics. As it 
is shown above the dual superconductor hypothesis predicts a particular phase diagram in the 
external electromagnetic fields. The numerical investigation of this diagram can be used for 
further checks of the dual superconductivity of the vacuum. Moreover, the value of the critical 
fields corresponding to the phase transitions (if exist) could be used to determine the closeness 
of parameters of the dual superconductor to the type-I/II boundary. 
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